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Molecular random walks in a fluid and an invariance group 
of the Bogolyubov generating functional equation 

Yuriy E. Kuzovlev 
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Abstract 



The problem of statistics of molecular random walks in a classical fluid is analyzed by means of the BBGKY 
hierarchy of equations reformulated in terms of the Bogolyubov evolution equation for generating functional of 
many-particle distribution functions. A proper equivalent set of correlation functions is introduced so that all they 

Sare integrable, vanish in statistical equilibrium, otherwise accumulate statistical information about history of collisions 
of a "molecular Brownian particle" (test molecule) with other molecules of the fluid. An exact evolution equation for 
generating functional of such correlation functions is derived. Then it is shown that time-dependent solution to this 
equation, as well as a properly defined generating functional of static thermodynamically equilibrium correlations, 
possesses invariance with respect to a definite group of transformations of independent variables of the functional, 
O . if density of the fluid (number of molecules per unit volume) is treated as one of the independent variables. Such 

invariance results in infinitely many exact relations between the correlation functions and probability distribution of 
path of the molecular Brownian particle. Even simplest of these relations suggest significant restrictions on a profile 
^. ■ of the path probability distribution, even without literal solving the BBGKY hierarchy. 

m 

(N 

Key words: diffusion, molecular random walks, Brownian motion, BBGKY hierarchy, Bogolyubov generating functional, 
^sj ' Bogolyubov functional equation, kinetic theory of fluids 

PACS: 05.20.Dd, 05.20.Jj, 05.40.Fb, 05.40.Jc 

o 

oo 
O 



1. Introduction 



Everybody knows that all thermodynamic and kinetic properties of a fluid are determined by thermal 
motion of its atoms or molecules. This motion by its nature is unbound and represents a "random walk" or, 
in other words, "Brownian motion". What kind of statistics characterizes it? 

Strangely enough, this principally interesting question never was a special issue of those fundamental 
approach to statistical mechanics of fluids what originated from the Bogolyubov's book [Tj and resulted in the 
Bogolyubov-Born-Green-Kirkwood-Yvon (BBGKY) infinite hierarchy of equations |1I2I3|4|5] , Unfortunately, 
next work of Bogolyubov and other theorists was concentrated at not conscientious solving these equations 
but forcible reducing them to a single "kinetic equation" for one-particle distribution function, at that being 
guided by the old Boltzmann's "molecular chaos" hypothesis. The latter assumes that two molecules entering 
into collision are statistically independent one on another, at least in case of dilute gas. 
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Such hypotheses rise up from naive identifying dynamical independence of particles and statistical one, 
in spite of that the first concerns individual phase trajectory of a dynamical system while the second 
statistical ensemble of trajectories corresponding to some ensemble of their initial conditions. But Krylov in 
his remarkable book [6] (first published in 1950) finally proved that dynamical independence gives no grounds 
for statistical one. Although, to tell the truth, this is clear as it is. Analogously, in a statistical ensemble of 
the mankind histories any two contemporaries would have somehow statistically dependent destinies even if 
they reside in different continents and do not influence one another. All the more if they do. Saying about 
gases, for example, in [Jj (see also [8]) it was shown that in spatially non-uniform statistical ensemble the pair 
distribution function for atoms in pre-collision states does not reduce to one-particle distribution function, 
even under the limit of infinitely dilute gas (or so-called Boltzmann-Grad limit). 

In fact, fortunately, the BBGKY hierarchy is self-consistent and has no need in extraneous hypotheses. It 
is quite another matter if a technique of its correct solving is in embryo. In respect to the above formulated 
question, an approximate approach was suggested in [7J (or see [8]). Its results were principally different 
from predictions of conventional theory operating with the Boltzmann and Boltzmann-Lorentz or similar 
equations 2 3 4 5J implied by the "molecular chaos" hypothesis. But since then no response to [7J was 
noticed. Hence, importance of the question still is not realized. This circumstance stimulated work [9J where 
due to simplification and improvement of method of [7] a whole probability distribution of Brownian path 
of a test atom was approximately found (instead of its four lowest statistical moments as in [7])- Recently in 
[TO] and then in [11112] a new approach to the same problem was developed which exploits the "generalized 
fluctuation-dissipation relations" |13|14j . A peculiar exact "virial expansion" for the probability distribution 
of Brownian path was obtained which allows to describe characteristic spatial profile of this distribution 
without any quantitative approximations at all. At that, results of [9] were confirmed. 

Well, clearly, now considerations of the BBGKY hierarchy are next in turn to be improved. With this 
purpose, in the present paper we make use of the Bogolyubov's formulation of the BBGKY hierarchy 
in terms of generating functional of many-particle distribution functions and evolution equation for this 
functional pQ. Notice that Bogolyubov himself analyzed stationary solution to this equation only, that is 
mere thermodynamically equilibrium functional, and in fact had not investigated its non-equilibrium time- 
dependent solutions. Apparently, because of this an interesting property of the generating functional and its 
evolution equation remained unnoticed, and we had to catch it. 

The matter concerns (i) logarithmic variational derivative of the equilibrium functional, (ii) generating 
functional of properly introduced non-equilibrium n-particle "correlation functions" which together repre- 
sent statistical correlations between a total previous path of the "molecular Brownian particle" (test atom) 
and current state of the fluid, and (iii) time evolution operator for this new functional. All these three 
objects turn out to be invariant with respect to definite continuous group of transformations of their inde- 
pendent variables. The latter are (i) arbitrary probe function defined in the 6D //-space and (ii) density of 
fluid (number of particles per unit volume) . 

It will be shown that this invariance property implies an infinite set of exact relations connecting the 
probability distribution of Brownian path and the correlation functions . A solution to the BBGKY hierarchy 
automatically satisfy all these relations. One of them coincides with the mentioned "virial expansion" . 
Therefore they present essential information about solution to the BBGKY hierarchy and even can prescribe 
it accurate to quantitative details. Let us go to these curious things. 

2. Equations of molecular Brownian motion and Bogolyubov's generating functional 

Let a macroscopic box with volume fi contains N 3> 1 identical atoms plus one more different but 
similar particle which will be under our special attention. The atoms have mass m , coordinates Tj and 
momenta Pj (j — 1,2... N) and all interact one with another by means of spherically symmetric potential 
<l> a (rj — rfc) = <f> a (\vj — Yk\/r a ) . The additional particle has mass M , coordinate R, momentum P and 
interacts with all atoms through spherically symmetric potential &b(Yj — R-) = 4>b(\^j — R-lAb) ■ Due to 
the interactions, this particle undergoes, along with atoms, a chaotic wandering, therefore we can name it 
"Brownian particle" (BP). We want to consider probability distribution of its position, R(t) , assuming that 



at some initial time moment, e.g. t = , it is certainly known: R(£ = 0) = Ro . 

On the contrary, personal coordinate of any concrete atom all the time is unknown and possesses nearly 
uniform probability distribution over tt . In such situation, normalized particular distribution functions for 
the BP together with any n = 0, 1, 2, ... atoms, following Bogolyubov [I], can be written as 

£>„(*, R,rWp,pW) = fi- n F n (t,R,rW,P,p(")|R ;^o) , 

where r(") = {ri... r„ } , pM = {pi... p n } , and new distribution functions F n depend on Q through mean 
density (number of atoms per unit volume), vq = N/Q , only. Argument Ro is added in order to remind of 
the mentioned initial condition for BP. At that 



F„(t,R,rW,P,p(")|R ;^o) = fi n / ... / Av(*,R,rW P,pW) 



n+l N 



with J, ... = J J ... drkdpk and Dn being normalized probability distribution function of the whole 
system. This distribution function (DF) is solution to the Liouville equation 
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with initial condition, in the framework of the Gibbs canonical statistical ensemble, 



D N (t = 0,...) = £(R-Ro) e 



-H N /T 



f rfR I dP f ... f S(R-R )e- HN/T 



N 



(1) 



Symbol H n hear means Hamiltonian of system " n atoms plus BP in the box" , L n + L„ ox ' is the Liouville 
operator corresponding to H n , and L n is its autonomous box- independent part: 



P d 

M "d~R 



E 



Pi 9 



E^+ E %> 



j=i 



l<j'<fc<n 



where 



L > X<1 ^ a -- i^-ap- 



L? fc = V* (rj - r fc ) • U- - -5— 



9pj 9p fc 



and V$(r) = d<P{r)/dr . 

Choosing initial distribution in the form ([I]) , we thus confine ourselves by investigation of random walk 
of the BP in a thcrmodynamically equilibrium fluid. 

Under reasonable interaction potentials (f> a ,b(x) (fast enough tending to infinity at x — * and to zero at 
x — > oo ) the particular DF F n have definite "thermodynamical limit" when walls of the box withdraw to 
infinity, N — > oo and £1 — > oo at iV/il = ;/o — const . In this limit the Liouville equation produces pQ 
infinite hierarchy of the BBGKY equations for F n 2 3 4 5j. In our case it appears as 



dF P dF , p 



dF\ 

dF n 

dt 



Lj. Fi + v / L\ F 2 +u L a 12 F 2 , 



= L n F n + vq I L n+1 F n+ i + vq y j \ L°[ n+1 F n+ i 

n+l i =1 n+l 



(2) 
(3) 
(4) 



The only difference of these equations from routine ones [1I2I3I4I5] is presence of a distinct particle, namely, 
the BP. Integration over its variables eliminates equation (2} and turns ([3])- ([4} to simplest BBGKY hierarchy 
for a system of equivalent particles. According to ([1]), initial conditions for our hierarchy are 

^o(0,R,P|R o ;^o) = S(R-Ro)Gm(P) , 

i^(0,R,r("\P,p( n >|R o ;^ o ) = S(R - R ) i^(r x ... r„|R; v Q ) G M (P) f[ G m (Pj) , 

where 

G m (p) = (27rTm)" 3/2 exp(-p 2 /2Tm) 

is Maxwell momentum distribution of a particle with mass m , and F„ (ri... r„|R; z/o) is usual thermody- 
namically equilibrium joint DF of coordinates of n atoms and BP corresponding to the canonical ensemble. 
Of course, eventually we are most interested in evolution of F (t, R, P|Ro; vo) and especially in probability 
distribution of BP'path, 



V {t,AR;is ) = f F {t,R,-p\Ro;v )dP (AR = R R 



It is important that F n by their definition naturally satisfy weakening of correlations between distant 
particles. Under initial conditions ([I]) and ([5]) this property concretizes to 

Fi(t,R 1 r 1 ,P,pi|R ;i4>) -» F (i,R,P|R ;i/ ) G m (pi) , 

^'(nlRo^o)- l , _ 

(6) 

F n (...r k .,.p k ...) -> F„_i(...r fc _i,rfe + i... Pk-i,Pk+t—)G m {pk) , 
F^(...r k ...) - F<?l(...T k - 1 , Tk+1 ...) , 

when ri — > oo and r^ — ► oo , respectively. 

It will be convenient, again following pQ, to collect all the DF into one generating functional 

OO 7) /* Z 1 

f{t,R,-p,Tp\R ;vo} = F (t,R,-p\Ro;vo) + J2^1 / - / ^(i,R,r (n) ,P,p (n) |R ;^) J] ^fo.p.,-) 
and correspondingly all the equations ©-(U) into single equation for this functional: 

dT t p dT ? ( s \ _ r Pl a &f , 



<9i M 9R \SipJ J m <9ri 5ip(xi) 
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(7) 



<ty(an) 2J J l ^ v ±/JL ^ v " yj i2 5il>(xi) 8il>(x2) 



For brevity here and below single argument aij deputizes for the pair {r^p.,} . 

Equation ([7]) becomes direct analogue of equation (7.9) from [1] if replace functional argument u(x) by 
i/oijj(x) and use obvious identity f t J 2 L\ 2 ... — [I]. Initial conditions ([5|) take the form 

^{0, R,P,^|R ;^o} = 5(R-Ro)GM(P)^" (e9) {0|R;^o} , (8) 

where generating functional for the equilibrium DF and new argument for this functional are introduced: 



f W W|R;"o} = 1 +£4 /- [Fi eq) (r 1 ...r n \R^ Q )f[cb(r 3 )dr 3 , 

71=1 ' ^ ^ j=l 

0(r) = / ^(r,p)G m (p)dp 



From physical point of view, expression ((5|) represents thermodynamical equilibrium, since fixation of start 
BP's position in no way disturbs it. Nevertheless, ([8]) is not a fixed point of equation |J7J|. The latter will be 
achieved at t — > oo only, when BP's position R(i) becomes completely uncertain. Hence, formal statistical 
equilibrium is presented by functional Gm(P) J-( eq ' {<j)\R; v$} , and we can write 
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G M (P)^ (e9) {</>|R;^} = o (9) 



This equation easy transforms into another one in terms of (/>(r) itself: 



d V$ 6 (r - R) 
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\J { l + m] V^-r)^^^ (10) 



This is analogue of equation (2.14) in pQ. As combined with boundary conditions at infinity ^ it quite 
unambiguously determines DF F n e and thus initial condition to equation ([7]). 

3. Correlation functions and historical correlations 

Though DF F n and thus F n (t = 0, ... ) include all equilibrium inter-particle correlations and at t > 
our system remains in equilibrium, nevertheless at t > some excess correlations do arise, so that 

F n (t,R > n...r B> P,pi...p n |R ;no) ^ F (i,R,P|R ; i/ ) F^(r 1 ...r n \R;v ) JjG m (pj) , 

3 

as opposed to what one could think. To describe a difference between left and right-hand side here, let us 
introduce "correlation functions" (CF) V n (t, R, ri ... r„, P,pi ... p„|Ro; i^o) and their generating functional 

V{t,R,P,V|Ro;^o} = Vb(t,R,P|Ro;^o) + X]4 /- /K(*,R,rW,P,pW|Ro;^) f[i>(r 3 , Pj ) 
and define all them in the language of functionals be relation 

F{t, R, P, ^ |R ; M = ^F {eq) {<f> |R; M V{t, R, P, </> |R ; M , (n) 

where </>(r) = f ip(r, p) G m (p) dp as in IjHJ). The term "correlation functions" is traditionally used in 
statistical mechanics for various corrections to equilibrium or quasi-equilibrium distributions (see e.g. [21314] ). 
By this definition, in particular, 

F (t,R,P|R ;^o) = y (t,R,P|Ro;^o) , 

F 1 (t,R,r 1 ,P,p 1 |Ro;^) = Fo(i,R,P|Ro;^o)F 1 (e9) (r 1 |R;^)G m (pi) + Vi(i,R,ri,P,pi|Ro;^o) (12) 
Initial conditions ([5]) transform to 

V n (t = 0,...) = ^,o<5(R-R ) , V{i = 0,R,P, ^|R ;M = <5(R-R ) , (13) 

while the boundary conditions ^ prescribe 

V n>0 (t,...r k ...) -> at r fc -» oo (14) 

It should be underlined that presence in V\ of start position of BP, Ro , is principally important. As we 
already mentioned, if Ro was unknown then then all the DF would be invariable, that is all V n> o would 
turn to zero. This remark highlights two peculiarities of excess correlations described by V n> o . First, they 
are essentially spatial correlations. Second, they connect current positions of atoms with not merely current 



position of BP but its total previous displacement, or path, R — Ro , during time interval (0, t) . By this 
reason we can name them "historical correlations" . Further we will see that their actual extent is very closely 
related to statistics of the BP's path R — Ro . 

Substitution of (jTTJ) to equation |0 yields, with accounting for identity ([9]), equation for the generating 
functional of correlation functions: 

£ + £-£-<(*£)v + '(*.*£)v. 

where action of operator C is described by expression |(7J) while action of new operator CJ by 

l+^( Xl )}[l + ij(x 2 )} L a 12 ; Y [ G m (p 2 ) 



(r 2 ) Sip(xi 

1 2 



So complicated expression is payment for simple initial condition Ijl3|> . After that, variational derivatives of 
(fT5JI at ip — give BBGKY equation in terms of the CF. The first of them looks quite similar to (|2]): 

dV P dV , f~ b 



at M OR ^J L ^ < 17 > 

l 

Since Vo = i*o , comparison of (jTTJ) with @ shows that quasi-equilibrium part of the pair DF Fi , i.e. first 
right-hand term in (JT3J), does not contribute to the collision term in @, and collisions of BP with atoms 
are wholly represented by pair CF V\ . The second BBGKY equation becomes 

^ = LiVi + P [Z b 2 V 2 + p [L a 12 [V 2 + F ( 1 eq \l)G m {l)V l {2)+F[ eq) {2)G m {2)V l {l)] + 

2 2 (18) 

+ G ro (l)|v$ 6 (R-r 1 )^ e9) (l)+ i , o y'v$ fe (R-r)a| e9) (r,r 1 |R^o)dr|^ + ^^ 

with F^ eq) (k) = Ft q) {r k \R^ Q ) , Vi(k) = Vi(t,R,r fc ,P,p fc |Ro;»b) and G m (k) = G m (p k ) . The function 

Cf^n^R; v ) = F 2 (e9) ( ri ,r 2 |R; v ) - F[ eq) (r^R v ) F[ eq \v 2 \R u ) 
is pair CF of the fluid under given BP's position. 

4. Equilibrium generating functional and its invariant transformations 

Next, consider more carefully the equilibrium generating functional T^ eq > . It is of interest for us here 
since operator (|16p includes its logarithmic variational derivative. Let us write the latter as 



— — - = v Q C{r,<f>\R; v \ 



VQ 



l + cf 9) (r|R;^)+ £ 5" / - / cS(r,r 1 ...r n |R;,„) 1]^)^ 

where C{ eq (r|R; v ) = F^' (r|R; i> ) — 1 , pair correlation function C 2 (r, ri|R; ^o) was presented just 
above, third-order correlation function or, to be precise, cumulant function is 

Ci eq \r,r u r 2 ) = F 3 (e<?) (r, n , r 2 ) + 2F 1 W(r)F 1 M(r 1 )F 1 M (r a )- 

_ F M (r)j pM( ri , r2 ) _ ^)( ri )F 2 M (r,r 2 ) - F^\t 2 )F^\t, Ti ) , 



for brevity omitting arguments R and is , and so on. Then equation (|10[) transforms to 

C<««>{r,0|R;i*} 



d V$ b (r - R) 



- / 1 + cj> r' V$ a (r' - r) x ,/ dr' + 

r./ <ty(r') ( 20 ) 



+ C^{r, |R; ^} Y j [l + 0(r ' } ] V$a(r ' " r) C(Gg) { r '' ^ l R ; **> dr ' 

This nonlinear equation can not be solved in a closed form but it will help us to perceive important features 
of functional C^ . 

With this purpose, first, notice that due to the DF's asymptotic at infinity © at any reasonable potentials 
of atom-atom and atom-BP interactions, similar to (1141). 



if r\ 



Ci e + 9 l(r,r 1 ...r n |R;^o) 

cf^OilRjz/o) -»■ if R-n->oo 



(21) 



-<( e ?) 



-l( e ?) 



C^^ri... r„|R; is {) ) - C^ 1 (r 1 ...r n |i/ ) if R - r, , - oo , 
where cumulant functions C„ (ri...r„|^o) describe fluid without BP and satisfy 



c£+{(r,ri...r n |i/ ) -»■ if r, - r -> oo 

Moreover, any of the cumulant functions C„> -^ri... r n |R; j/q) tend to zero at Yj — ► oo in so fast way that 
all they are integrable in respect to ri... r„ . 

Consequently, functional C^ eq ^ keeps a sense when <fi(v) turns into nonzero constant, <j>(r) -^ a = const , 
as well as when one replaces (f>(r) by a + (f>(r) with a = const . Thus we can introduce function 



C(a,is ) = lim C^{r,a\R;u } = 1 + V ^- / ... / c£g(r,n...r n |m) dv x ...dv n , 

R r -^oo ^ — ' Til I I 

n=l J J 

where, clearly, integrals factually do not depend on r . Besides, introduce new functional 

C^Uv ,61 m = C (eg) {r,a + 0|R;^o} C^>{r,q + 0|R; ^ } 

ff t '*' ' limR.r^CMjr.alRji/o} C(a,is ) 

Under change 4>{r) — > a + 0(r) equation (f2TJ)) evidently turns to equation for C^ {r, |R} as follows, 

. V$Q(r r ^[0(r')/(l + a)] dr 
V$ (r'-r)C( e 9){r',^|R}dr' 



(22) 



(23) 
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+ C^){r,0|R} 



C^){r,0|R} = I J 
u C(a,v ) (1 + cr) 
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1 + ff 



(24) 



Second, a simple analysis of structure of equation (|20j) shows that it uniquely determines its solution in 
the form of expansion (|19p if coefficients of the expansion satisfy conditions (|2"Tj) . 

Third, coefficients of similar series expansion of the functional (|23|) satisfy absolutely similar conditions, 
because, in consequence of ([21~j) . 
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(eg) 

n+/c+ 



1 (r, ri...r„, pi...p fe |R;i/ ) dpi— dp* 



if r 9 



4 e9) {r,0 = O|R} -» 1 if R-r 



(25) 



Fourth, the only formal difference of equation (j2"4")) from (j2"0|) is mere change of its two independent 
variables, concretely, <^>(r) — ► </>(r)/(l + a) and r/ — » is (a, isq) , where 



z/(ct, Vq) = is C{a,is Q )(l+a) 



(26) 



These observations taken together mean that a proper solution to equation ([M)l . C& {r, <j> |R} , is nothing 
but C^ eq ^{r, 4>/(l + <t) |R; v(<j, vq)} , that is, in view of (j23f and (p6|) . for any a = const identity 



1 + cr [ 1 + (T J 



^^{r^ + ^R;^} = -A^ C (^) r, -^— | R; K^o) (27) 

takes place. In essence two latter formulas define a group of such transformations of equilibrium generating 
functional C^ eq '{r, (f) | R; i>q\ which do not change its value: 

f( CT )C^){r^|R;^ }= ^ ^ C (eg) (r, ^ - 1 [ R; u(a, v Q ) \ = C^){r^|R; u } , 

(l + a)v Q I 1 + ct J 

f (a 2 ) f ( CTl ) - f (a! + <t 2 + axa 2 ) , ( 28 ) 

z/(ct 2 , ^(cti,i/ )) = v(o x +(T 2 + cricr 2 ,^o) 

At that, obviously, the group parameter should satisfy a > — 1 , and similar restriction should be imposed 
upon 4>(r) . But this restriction disappears if write a = exp (a) — 1 , so that in terms of new parameter 
T(a 2 )T(ai) = T(ai+a 2 ). 

Infinitesimal form of functional identities (I27|) or (|28|) looks better in terms of particular cumulants: 



(29) 



\uk(v) + [l + x(i/)]i/^}[C^(ri...r n |R;i/)+* n)1 ] = v J C^{(n...r n> r'|R;i/) dv 1 , 
= v [ C { 2 eq) (r,Q\v)dr 

CT=0 J 
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'dC{a,v)~ 
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Recall that 




1 + x{y) 


f du 
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with P being pressure, is isothermal compressibility of a fluid [314115] . 

5. Invariance group of the historical correlations functional and virial relations 

Now, turn our attention to evolution equation (|15[) for the generating functional of correlation functions, 
V{£, R, P, "0 |Ro; vq} j an d consider its principal properties implied by structure of operator £ + £' , above 
mentioned properties of factor Jln^^l^R; vo}/5(j)(r) = ^C' e9 '{r, 0|R; vq\ which enters £' , initial 
condition (fT3|) and, of course, asymptotic boundary conditions (fT4]) . 

What is important, the initial condition for V does not include independent variables ip = ip(r,p) an d 
ua at all. Therefore all effects of their transformations are purely determined by properties of £ + £' 
itself under boundary conditions (|14[) . One role of the latter is as just above: they allow to extend functional 
V{t, R, P, ip |Ro; vq} to argument a + ip(r, p) with a = const in place of ^(r, p) . Another role of conditions 
(fT4"|) is that due to them, and according to definition ([7]) of operator £ , variable ip(xi) m the first term of 
£ V can be shifted by arbitrary constant: 



fit \ Pl ° SV [\ j. U M Pl d SV 
/ V(X\) tt- 777 — 7 = / [v + V{xi)\ 77— T-77 — 7 , 

J m OTi oyj(xi) J m OYi Olp(Xi) 



Sip(xi) J m 9ri 5ip(x\) 

1 1 

where a — const , e.g. a — 1 . This is principal difference of expression £ V from £ T . As the consequence, 
taking in mind action of £ onto V , one can write 



Visual study of operator L' defined by (|16p. with taking into account identity (|2"7)l for i/oC' e, '{r, </>|R; fo} , 
shows that C possesses the same invariance property, if transformation of i/)(r,p) is accompanied by 
transformation of the density argument vq in conformity with (I26|) and (j22[) : 



z \""' + *-h) " '''(^'"^' wm ) (31) 

Because of pop and pip equation Q15p under initial condition Q13p implies similar invariance property for 
its solution: 

V{t,R,P, a + ^|R ;r/ } = V It, R,P, J^- | Rq; i/(o-,*b)j (32) 

with arbitrary cr(r) = const . In other words, the generating functional V of time-dependent "historical 
correlations" also is invariant under above mentioned group of transformations: 

f (ct) V{t, R, P, V I Ro; H>} = V{t, R, P, \^- - 1 1 R ; v(<r, is )} = V{t, R, P, V I Ro; H)} , (33) 

where left identity, as combined with (|2l)|) and (|2"2"j) . defines action of the group onto V . 

This is main formal result of the present paper. It would be difficult to obtain this result from a chain of 
equations for particular DF, ©-((I]) or all the more (|T7|) . (|18p . etc. Thus, perhaps, for the first time a real 
profit of Bogolyubov's generating functional for non-equilibrium DF has been demonstrated. 

Infinitesimal form of (l33l reads 



|nx(i/) + [l + x(i/)]i/^-| V n (t,R,rW,P )P W|Ro;i/) = 

= v J K+i(t,R,r(" +1 \P,p( n+1 )|R ; v) , 

n+l 

thus at n > connecting derivative of any BP-atoms correlation in respect to density of atoms and next- 
order correlation integrated over variables of additional outer atom. At n = we connect derivative of the 
very BP's probability distribution and integrated pair historical correlation: 

^ ^o(t,R,P|Ro ;i /) /■„,.„ p IR v 

[l + M 1 ')] ^ = / Vi(i,R,ri,P,pi|Ro; !/) (35) 

1 

with x{v) defined in ([29]) . 

From the other hand, if considering ^(r, p) as infinitesimal parameter instead of a and performing in p2p 
fc-order variational differentiation over ip at point ip = , one comes to rather different type of relations. 
At k > , 



y fc (i,R,r( fc ),P,p( fc )|R ;K^^o)) = 

(1 + crj^oj 

00 y n r7 n r r (36) 

= %(i,R,r«,P,p( fe )|R ;^) + E^T- / - / ^•+n(*,R,r (fe+ " ) ,P, P (fc+ " ) |Ro; *o) 

™ _i fc+1 fc+ri 

For fc = , i.e. merely under putting on tp = , identity (|32l) yields 



Vb(i,R,P|R ; vfovo)) = 7 (i,R,P|Ro; ^0) + ^ ^T 7 "7 t r »(*> R > r(n) . p >P ( '°l B <>; »*) (37) 

71 = 1 1 n 

Of course, all terms in this formula actually depend on the BP's path AR = R — Ro only. 

Formulas (|37|) and ([36]) . by their sense, can be interpreted as a sort of virial expansion of the BP's 
probability distribution and BP's correlations with fluid, with those peculiarity that in opposite to usual 



virial expansions of thermodynamic quantities [15j or kinetic coefficients |16] our ones are expansions over 
relative change of density, v/vq — 1 , rather than its absolute value. Then corresponding differential forms 
(l35l) and (l34l) can be named "virial relations" . 



6. Some principal consequences of the virial relations 

Although invariance property (|32|) of the functional V does not substitute solving the equation (jf 5p 
of its evolution, corresponding virial relations bring important information about structure of correlations 
functions and, eventually, character of evolution of the BP's probability distribution 

V (t,£R,P\v) = Vb(t,R,P|Ro; v) 



For example, let us consider relation (|35[) supposing for simplicity that fluid density v in ([35)1 corresponds 
to a dilute gas: 47rr^/3 <C 1, 47rr^/3 <C 1 (characteristic radii of interactions r a and r b were introduced 
above). Then x{v) <C 1 and, as it is known (see e.g. [16]), mean free path of our molecular-size BP and mean 
free path of atoms both are inversely proportional to gas density: Ab ~ (nr1v)~ l , A a ~ (^r^v)" 1 , together 
with diffusivity of BP, D b = D b (v) ~ u&Ab oc v^ 1 , and diffusivity of atoms, D a (v) ~ v a A a oc v~ x , where 
Vb = \JtJM and v a = y/T/m are thermal velocities. Besides, we take in mind large enough temporal and 
spatial scales: |AR| ^> Af, and t ^> 77, with r& ~ A^/vf, being mean free-flight time of BP. 

First, it is useful to compare relation ([55]) . along with (154")) . and relations (|29[) . In the latter, in general, 
max|Cn I ~ 1 while characteristic scale of equilibrium inter-atom correlations equals to r a . Therefore 
right-hand sides of (|29|) can be roughly but surely estimated as 

Jc^\(T l ...r n ,r'\R;u)dv' ~ r 3 C^>(i 1 ...i n \R i u) 
On the left, consequently, 

^^(r^r^R;,) ^ ^ a^)( ri ...r n |R;^ « <#«>(,!... r n |R;„) , 
in particular, x(i') <C 1 . In opposite, in (|34 | - ([35|) we can at first roughly but surely estimate left sides: 

i , av.(.,B 1 ^P,p<->|B.,» „ K>(t]RirWiP]pW|B<ji „ ) 

Indeed, a relatively small variation of gas density results in similar relative variation of BP's diffusivity and 
therefore in comparable relative change of probability of its displacement, as well as connected probabilities 
and correlations. For instance, if the hypothetical Gaussian asymptotic [3I4I5J of BP's distribution at t 3> r b 
and |AR| ~ %/£>&£ » At realized, 

Vo(t,AR,P|i/) -» Vfc(i,AR,P|i/) = [47rD 6 (j/)t]- 3 / 2 exp hAR 2 /4Z? 6 (^)i] G M (P) , (38) 

then one would have, evidently, 

3 AR 2 



— w Vo(t,AR,P|i/j 



tf 



;/ 



2 AD b (v)t 



(39) 



Hence, for right-hand sides of (f34|) and ([35]) we obtain estimates J n+1 I4+i ~ dV n /dv ~ 1/ 1 I/ n and 



I^^Rr^Rp^Ro; „) ~ ^o(f,AR,PM ^ ^ y^^^ 



(40) 



with coefficient i/ - 1 which is (r^) - : 3> 1 times greater than in case of equilibrium correlations. 

This observation prompts that characteristic volume, Q c , occupied by a historical correlation between 
BP and an atom by order of magnitude is equal to v~ 1 , that is volume displayed per one atom of fluid. 
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Physically, such statement follows also from almost trivial syllogism: 

(a) historical correlation between BP and atom springs from their current (or soon forthcoming or just 
happening) collision; 

(b) collision is (or will be or was) possible among those particles only whose relative position vector ri — R 
belongs to the "collision cylinder" ; the latter has radius rs r b , is oriented along relative velocity of BP and 
atom, Vi — Vo = pi /m — Po/M , and its length in this direction has an order of A;, (assuming for simplicity 
that r a ~ r b and A a ~ Ab ) since collision of stronger separated particles is prevented by other particles; 



(c) consequently, volume of the pair correlation fl c 



r 2 b A b 



Formally, such statement, i.e. fi c 



max |Vi(t,R, ri,P,pi|R ; v)\ 

ri 



follows from (j40|) if one presumes that 



dy (i,AR,Pk) 



ch 



G m (pi 



y (t,AR,P|i/)G? m (pi) 



where maximum is achieved somewhere at |ri — R| ~ r b . This presumption seems likely from the point of 
view of equation (fP8)) and well agrees with considerations of pair correlations in conventional theory [31415] . 
However, the resulting conclusion that fl c ~ v~ 1 contradicts conventional models of molecular diffusion 
("hybrid" models, in terms of 0], as they combine dynamics and stochastics) . In order to see this, let us 
consider relation (fl"2"l) (in fact, definition of Vi ) taking |ri — R| > r b . Since left side of (JT3J) represents a 
probability density, it is certainly non-negative. Therefore everywhere in the mentioned region, including 
most part of the collision cylinder, 



V 1 (t,R,ri,P,Pi|Bo;i') > -V (t, AR,P| v) G m (p x ) 
Hence, according to the estimate il c ~ v~ 1 , 



/Vi(t,R,r 1 ,P,p 1 |R ;i/) > - n c V (t, AR,P| v) ~ - v~ x V (t, AR,P| v) 



(41) 



(42) 



regardless of t and AR . This inequality as combined with (f3"5"|) or (|40[) and confronted with (|3"5)) clearly 
shows that formula (|38)) can not represent a true asymptotic of the BP's probability distribution! 

A more rigorous substantiation of this striking statement was suggested in [11112) . There (see also |10j ) 
it was shown that true asymptotic of Vo(t, AR, P| v) as a function of AR possesses power-law long tails 
cut off at |AR| ~ Vbt . This qualitative prediction agrees with quantitative result of approximate analysis 
of the BBGKY hierarchy undertaken in [S] : 



V (t,AR,P\u) 



T(7/2) 



[4nD b (v)t 



3/2 



1 



AR 2 



4,D b (v) t 



■7/2 



e 



|AR| 
v b t 



Gm(P) , 



(43) 



where function 0(a;) w 1 at x ^C 1 and in a fast way tends to zero at x > 1 . An origin of such statistics 
of equilibrium molecular Brownian motion was discussed more than once as long ago as in [17I18J and 
later in [71819119] (see also references therein). In short, such statistics says merely that dynamics of a very 
many-particle system always is unique and can not be imitated by primitive stochastic processes. 

7. Conclusion 

To resume, first, in this paper we rigorously formulated the problem of random walk of small molecular- 
size Brownian particle (BP) in classical thermodynamically equilibrium fluid. We considered correspond- 
ing BBGKY hierarchy in terms of Bogolyubov type generating functionals and derived a linear evolution 
equation for generating functional of time-dependent "historical correlations" which together accumulate 
statistics of collisions of BP with fluid atoms and eventually determine statistics of BP's path. 

Second, we showed that both the evolution equation and its solution, along with generating functional 
of static equilibrium correlations, are invariant with respect to definite continuous group of transformations 
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of their independent variables including fluid density. As the consequence, we found an infinite set of orig- 
inal exact "virial relations" which connect the BP's path probability distributions and various correlation 
functions at different values of the density. 

Notice that particular relation ([37]) first was found in [10] in the form integrated over BP's momentum 
and in [ll|12j in the full form. There it was derived from most general properties of the Liouville operator, 
as an example of the "generalized fluctuation-dissipation relations" . 13:14] ■ Here, formula ([37]) appeared, 
in company with infinitely many new relations (|34[) and (|36[) . as consequence of hidden symmetry of the 
BBGKY hierarchy. In principle, of course, these are allied approaches. Nevertheless, they use very different 
formal techniques, therefore confirmation of previous results certainly is useful. The more so as it produces 
much wider results which give a new sight of the BBGKY theory. 

Third, we demonstrated that virial relations bring significant information about solutions to the BBGKY 
hierarchy even without literal solving it. In particular, they impose principal restrictions on possible asymp- 
totic profile of the BP's path probability distribution, definitely forbidding the Gaussian profile implied by 
invented stochastic or hybrid models of random walks [31415] , in agreement with results previously obtained 
by approximate methods in [7] (or see [8]) and in [9]. 

Undoubtedly, this is not all of the benefits from the invariance group of BBGKY hierarchy represented by 
virial relations. It can extended, from one hand, to non-equilibrium random walks at presence of an external 
force acting onto BP. From the other hand, to more usual problems when probability distributions under 
interest describe not information on position of some distinct particle but hydrodynamic fields of initially 
non-equilibrium or externally disturbed fluid. Curious feature of the virial relations (|36|) - (137|) is that they 
establish contacts between states of a fluid under arbitrary different values of pressure, e.g. between dense 
liquid and dilute gas. May be this will help in constructing approximate solutions to the BBGKY equations. 
Anyway there is a lot of interesting tasks for the future. 
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